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SECTION  I 


INTRODUCTION 

The  problems  of  testing  the  hypotheses  on  the  structures  of  the 
covariance  matrix  of  the  multivariate  normal  population  have  received 
considerable  attention  in  the  literature  since  these  problems  have 
applications  in  various  disciplines.  Wilks  [1,2]  is  one  of  the  earliest 
workers  in  this  area.  In  this  report,  we  investigate  approximations  to 
the  distributions  of  the  likelihood  ratio  statistics  for  testing  the 
hypotheses  of  (1)  multiple  independence  of  several  sets  of  variables, 

(2)  sphericity,  (3)  equality  of  the  covariance  matrix  to  a given  matrix, 
and  (4)  equality  of  the  covariance  matrices  of  independent  sets  of 
variables.  These  approximations  are  based  upon  fitting  suitable  Pearson 
type  distributions  by  using  the  first  four  moments  of  the  above  statis- 
tics. Using  the  above  approximations,  we  construct  percentage  points  of 
these  statistics.  The  accuracy  of  these  approximations  is  sufficient 
for  practical  purposes. 

Exact  distributions  of  the  test  statistics  considered  here  are  very 
difficult  to  compute.  One  can,  of  course,  use  Box's  asymptotic  ex- 
pression by  taking  a sufficient  number  of  terms  to  compute  the  tables 
and  this  may  be  preferable  to  exact  expressions  from  a computational 
point  of  view.  However,  this  asymptotic  expression  is  also  difficult 
to  compute  when  we  need  to  take  several  terms  in  the  series  and  this  is 
the  case  when  vhe  sample  site  is  not  large.  Since  Pearson  type  approxi- 
mations are  quite  simple  from  a computational  point  of  view  and  since 
their  accuracy  is  sufficient  for  practical  purposes,  they  are  definitely 
preferable  to  either  exact  expressions  or  Box's  asymptotic  series  if  we 
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are  interested  in  computing  the  percentage  points  of  the  statistics 
considered  in  this  paper. 
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SECTION  II 


STATEMENT  OF  THE  PROBLEMS  AND  PRELIMINARIES 


Let  X1  * (X' , . • . ,X')  be  distributed  as  a multivariate  normal  with 
- -i 

mean  vector  y'  and  covariance  matrix  I.  Also,  let  E(X^)  * y and 

q 

e{ (X^  - y^)  (Xj  - y ^ ) * } * E^,  where  X^^  is  of  order  pA  x 1 and  s » Z p^. 

In  this  paper,  we  consider  approximations  to  the  distributions  of  the 
likelihood  ratio  statistics  for  testing  the  hypotheses  H , H4# 


and  Hg  where 


«1!  “ 0 (i  + j - i#...#q) 


H.  : E - o2  E (o2  is  unknown,  Z is  known) 
2 o o 


H Z * E 
3 o 


H.  : Z. . - ...  » Z (under  the  assumption  that  H,  is  true 
4 11  qq  1 


and  p,  * ...  ■ p ) 
1 <J 


E * ...  - E 

u V*1 


«5  ’ 


EVl'Vl  “ 

Z * , . » . , * ...  * I 

‘S-l*1'^-!*1  q'q 


j 


Here 


q*  « o,  q*  * Z q±  and  q*  - q. 


i-1 


These  distributions  are  approximated  with  Pearson  type  distributions. 
A brief  description  of  the  family  of  Pearson  type  distributions  is 
given  below. 
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If  £ (x) , (c  < x < d) i is  the  probability  density  function  be- 
longing to  the  family  of  Pearson  type  distributions,  then  f(x)  satisfies 
the  differential  equation 

(bQ  + bx*  + b2x2)  ^LL  - (x  + a)f  (x).  (1) 

Let  m£  ■ (x)dx  and  m^  - (x  - m^)*! (x)ax.  It  is  known  that 

2 

a - m3(»4  + 3m2)/A 
b0  “ ■b2<4“2"4  ' 

bj  - -m3<«4  + 3m2) /A  (2) 

b2  “ *(2n2m4  ' 3*3  ' 6m2)/A 

A ■ 10mm.  - 18m2  - 12m2 
2 4 2 3 

The  family  of  Pearson  type  distributions  is  completely  specified  by  its 
first  four  moments  and  the  type  of  the  distribution  can  be  determined  by 
examining  the  quantity  k where 

B1(B2  + 3)2 

* " 4(28j  - 3B1  - 6)(4B2  - 3BX)  ' <3> 

2 3 2 

0,  ■ m and  0.  - m /■» . If  ic  is  negative,  we  classify  it  as  tne 
1 3 2 2 4 2 

Pearson's  Type  I distribution.  When  ic  is  positive,  it  is  classified  as 
the  Pearson's  Type  IV  or  Type  VI  accordingly  as  ic  is  less  than  or  greater 
than  unity,  respectively.  For  a discussion  of  Pearson  type  distribu- 
tions, see  [3].  Tables  for  the  percentage  points  of  Pearson  type  dis- 
tributions are  given  in  Johnson,  Nixon  and  Amos  (4)  for  some  values  of 
the  parameters. 
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SECTION  III 

MULTIPLE  INDEPENDENCE  OF  SEVERAL 
SETS  OF  VARIABLES 


Let  « (X^,...,X^),(j«l,...,N),  be  jth  independent  observation 
on  X'  and  let 


where 


A, , • • • A, 
11  lq 


A _ _ . . . A_ 
21  2q 


A ....  A 

qi  qq 


N 

A*“  ' V’ 


and  X^#  • — I Xtj.  Then,  the  likelihood  ratio  statistic  for  testing 
is  known  to  be 


1 q . 

n |a 

j-i 


11 ' 


(4) 


The  nonents  of  the  statistic  are  given  by 


!<vj) 


n r fi(n+i-i>  + hi 

q 

n 

nArf  i(n+x-j>  1 

1 

i-i  *-2  J 

i-l 

l i-i  L 2 J 

t 

1 * 

Pi 

n 

11 

n r 

±-<n+l-j>  ♦ h 

J i-l 

j-i 

L 2 

(5) 


S 


where  n * N-l  and  T(n)  is  the  complete  gamma  function.  The  statis- 
tic V and  its  moments  were  derived  by  Wilks  [2] . According  to  the 

* < 

likelihood  ratio  test,  v*  accept  or  reject  accordingly  as  > c, 

* 

where  * -2 log  and 

r[v*  - (1  - a)  (6) 


The  distribution  of  is  quite  skew.  Hence,  we  approximate  the  distri- 
bution of  vj*  with  a suitable  Pearson  type  distribution,  where  b is  a 
properly  chosen  integer.  The  type  of  Pearson  distribution  that  has  to 
be  fitted  is  determined  by  computing  K given  by  Eq.  (3)  after  replacing 

ml'  m2'  m3'  m4  with  ^bl'  ^b2'  ^b3  ^b4'  resPectively»  ^ere  y£h 

is  the  h-th  moment  of  By  the  above  method,  it  is  found  that  the 

Pearson's  Type  I curve  fits  the  distribution  of  The  density  of 

the  Pearson's  Type  I distribution  is  given  by 


f(x) 


B(a+l,  e+lHa-a  ) 
o 


-<x-cr  )a(a-x)e 

a+e+l  o 

a < x < a 
o — — 


(7) 


where  a , a,  0 and  e depend  upon  a°,  b°,  b°  and  b°.  The  parameters 
o o 1 2 

o o o o 

a , b , b and  b are  respectively  equivalent  to  a,  b , b_ , and  b. 

0 12  0 12 

after  replacing  m£  with  In  this  report  we  consider  the  case  where 


Box  [5]  derived  asymptotic  expressions  for  a class  of  likelihood 
ratio  test  statistics  in  multivariate  statistical  analysis.  This  class 
includes  the  statistics  associated  with  testing  the  hypotheses  H^,  H2# 
H^,  and  H^.  The  number  of  terms  given  by  Box  is  not  sufficient  tc  get 
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the  desired  degree  of  accuracy  in  several  practical  situations.  So,  the 
authors  gave  terms  up  to  0(n~15)  in  the  appendix  applying  the  method  of 
Box. 


Consul  [6]  and  Mathai  and  Rathie  [7]  derived  exact  expressions  for 
the  distribution  of  V^,  and  these  expressions  are  very  diffucult  to 
coelute.  So,  the  authors  studied  approximations  to  these  distributions 
with  suitable  Pearson  type  distributions. 

Table  I gives  a comparison  of  values  obtained  by  the  Pearson  type 
approximation  and  the  asymptotic  expression  of  order  n 13 . In  Table  I, 
is  the  value  of  a if  we  use  the  Pearson  type  approximation  whereas 
is  the  value  of  a when  we  use  the  asymptotic  expression  of  order 
n 13 , wher*  a is  given  by 

plvi-cJHiJ  “ <1_a)-  (a) 

Davis  and  Field  [8]  computed  the  percentage  points  of  -2p  log 
for  some  values  of  the  parameters  by  using  the  Comish-Fisher  type  in- 
version (see  Davis  [9])  of  Box's  asymptotic  series  when 


In  Table  II,  the  entries  under  the  columns  L-C-K  are  the  values  of  c^ 
obtained  by  the  authors  with  the  Pearson  type  approximation  whereas  the 
entries  tinder  the  columns  D-F  are  the  corresponding  values  obtained  by 
Davis  and  Field  [8] . 
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Tables  I and  II  indicate  that  the  accuracy  of  the  Pearson  type 
approximation  to  the  distribution  of  is  sufficient  for  practical 
purposes.  Hence,  using  the  Pearson  type  approximation,  we  computed  the 
values  of  c^  for  p * 1,  2,  3;  q = 3,  4,  5?  a = 0.01,  0.05,  0.10;  and  M ■ 
1(1)20(2)30,  where  M = n - s - 3.  These  values  are  given  in  Table  III. 
When  q * 2,  it  is  found  that  the  Pearson  type  approximation  is  quite 
satisfactory,  and  the  results  are  reported  in  a companion  paper  [10] 
by  the  authors. 


TABLE  I 

COMPARISON  OF  THE  PEARSON  TYPE  APPROXIMATION 

WITH  THE  ASYMPTOTIC  EXPANSION  OF  ORDER  n'13 
FOR  THE  DISTRIBUTION  OF  vj 


U 

EM 

q=5 

P-1 

n 

ci 

ai 

°i 

ai 

a2 

10 

1.913 

0.05 

0.0499 

4.978 

0.05 

0.0488 

10 

2.778 

0.01 

0.0100 

6.333 

0.01 

0.0095 

15 

1.187 

0.05 

0.0500 

2.947 

0.05 

0.0497 

15 

1.723 

0.01 

0.0100 

3.742 

0.01 

0.0099 

20 

0.860 

0.05 

0.0500 

2.099 

0.05 

0.0499 

20 

1.249 

C.Oi 

0.0100 

2.663 

0.01 

0.0100 

30 

0.555 

0.05 

0.0500 

1.333 

0.05 

0.0500 

30 

0.805 

0.01 

0.0100 

1.690 

0.01 

0.0100 
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entries  in  this  table  are  the  values  of  c,  where 


s 

l 

: 

f 
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SECTION  IV 


LIKELIHOOD  RATIO  TEST  STATISTIC  FOR  SPHERICITY 


The  likelihood  ratio  statistic  for  testing  H.  : E = o E , where 

2 o 

2 

a is  unknown  and  is  known,  is  given  by 


(9) 


where  A was  defined  in  Section  III  and  tr  B denotes  the  trace  of  B.  The 
h-th  moment  of  is  given  by 


E(vh)  = 

2 r<sh  + 


s 

n 

i=l 


r [j(n+l-i)+h] 
r[j(n+l-i)  j 


(10) 


The  statistic  and  its  moments  were  derived  by  Mauchly  [11].  Using 
the  first  four  moments,  we  approximated  the  distribution  of  with 

the  Pearson  Type  I distribution,  where  b is  a suitably  chosen  integer. 
Expressions  for  the  exact  distribution  of  were  given  by  Consul  [6] , 
Mathai  and  Rathie  [12],  and  Nagarsenker  and  Pillai  [13],  but  these 
expressions  are  very  difficult  to  compute.  In  Table  IV  the  values  under 
the  column  LCK  are  the  values  obtained  by  the  authors  using  the  Pearson 
type  approximation  whereas  the  corresponding  exact  values  are  taken  from 
Nagarsenker  and  Pillai  [13] . This  table  indicates  that  the  accuracy  of 
the  Pearson  type  approximation  is  sufficient  for  practical  purposes. 

The  values  of  a in  this  table  are  given  by  the  relation 


P [V2^2IH2] 


(1  - a) 
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TABLE  IV 


COMPARISON  OF  THE  PEARSON  TYPE  APPROXIMATION  WITH 
EXACT  EXPRESSION  FOR  THE  DISTRIBUTION  OF  V2 


s=4 

s*5 

3=7 

n 

a 

L-C-K 

Exact 

L-C-K 

Exact 

L-C-K 

Exact 

6 

0.05 

0.0169 

0.0169 

0.0013 

0.0013 

— 

— 

6 

0.01 

0.0050 

0.0050 

0.00022 

0.00022 

— 

— 

iO 

0.05 

0.1297 

0.1297 

0.0492 

0.0492 

0.0029 

0.0030 

10 

0.01 

0.0726 

0.0726 

0.0242 

0.0242 

0.0010 

0.0010 

15 

0.05 

0.2812 

0.2812 

0.1608 

0.1608 

0.0368 

0.0368 

15 

0.01 

0.1967 

0.1966 

0.1052 

0.1052 

0.0207 

0.0207 

21 

0.05 

0.4173 

0.4173 

0.2877 

0.2876 

0.1111 

0.1111 

21 

0.01 

0.3264 

0.3264 

0.2156 

0.2156 

0.0761 

0.0761 

33 

0.05 

0.5833 

0.5833 

0.4663 

0.4663 

0.2665 

0.2665 

33 

0.01 

0.5013 

0.5013 

0.3910 

0.3910 

0.2125 

0.2125 

41 

0.05 

0.6507 

0.6508 

0.5453 

0.5453 

0.3515 

0.3515 

41 

0.01 

0.5769 

0.5769 

0.4741 

0.4741 

0.2939 

0.2939 
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SECTION  V 


LIKELIHOOD  RATIO  STATISTIC  FOR  TESTING  THE 
HYPOTHESIS  THAT  THE  COVARIANCE  MATRIX 
IS  EQUAL  TO  A SPECIFIED  MATRIX 

Consider  the  probler*  of  testing  the  hypothesis  H : Z ~ Z where 

3 0 

is  specified.  The  likelihood  ratio  statistic  for  testing  and  the 
moments  of  this  statistic  were  derived  by  Anderson  [14] . The  modified 
likelihood  ratio  test  statistic  (obtained  by  changing  N to  n in  the 
likelihood  ratio  statistic)  and  its  moments  are  as  given  below: 


V3  = (e/n) 


sn/2 


A Z 


-1 


n/2 


etr 


4 A 


(ID 


E(V^)  = (2e/n)Shn/2  \Z  |nh/2 


|l  + h Z 


I -n(l+h)/2 


II  r £y(n+nh+l-i)J 

n rj^j(n+i-i) 


(12) 


Using  the  first  four  moments,  we  approximated  the  distribution  of  with 

the  Pearson  Type  I distribution  as  in  Section  III,  where  b is  a suitably 

chosen  integer.  Korin  (15)  obtained  an  asymptotic  expression  of  order 
-15  * 

n for  the  distribution  of  * -2  log  Vy  Using  this  expression,  he 

★ 

computed  percentage  points  of  for  some  values  of  the  parameters. 

Nagarsenker  and  Pillai  [16]  obtained  an  expression  for  the  dis- 
* 

trihution  of  V^,  but  this  expression  is  complicated  from  a computational 

point  of  view.  Using  this  expression,  they  computed  exact  percentage 

* 

points  of  . In  Table  V we  compare  our  values  (given  in  the  column 
L-C-K)  obtained  by  using  the  Pearson  type  approximation  with  exact 
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values  of  Nagarsenker  and  Pillai,  and  the  values  obtained  by  Korin  [15] . 
Table  V indicates  that  the  accuracy  of  the  Pearson  type  approximation  is 
sufficient  for  practical  purposes. 
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COMPARISON  OF  THE  PEARSON  TYPE  APPROXIMATION 
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SECTION  VI 


LIKELIHOOD  RATIO  STATISTIC  FOR  TESTING 
THE  HOMOGENEITY  OF  COVARIANCE  MATRICES 


Let  X'  = (X X')  be  as  defined  in  Section  II.  Let  X. . 

-1  -q  -13 

(j=l,...,Nj  be  j-th  independent  observation  on  X^.  Also,  let 
Pi  = P(i  ■ 1# • • »q>  and 


N. 

l 

A.  . = Z (X. . - X.) (X. . - X.) 
3-1  jssl  -*3  -1  -1 


Wilks  [1]  derived  the  likelihood  ratio  statistic  for  H,  : Z,,  = ...  = Z 

4 11  qq 

when  Z = 0(i^j-l,..q)  and  derived  its  moments.  Let  n^  = N^-l,  and 

q 

n = Z n..  The  modified  likelihood  ratio  statistic  V (obtained  by  inter- 
i=l  X 

changing  wi*.h  n^  in  the  likelihood  ratio  statistic)  for  and  the 
moments  of  V are  given  below: 


V4  = 


q 

n /2 

n 

A 

g 

s ei 

gg 

nPn/2 

g=l 


gg 


n/2  q pn  /2 
n n 9 
9=1  9 


E(V‘,‘>  = I nphn/2/  n n 
/ 9=1 


/ 


q phn  /2 


(13) 


P 

“ 

[ q 

• n 

1 R r ( (n  +hn  +l-i)/2)/n  (n  +l-i)/2] 

i-1 

[g=i  9 9 9 

• | f l ( (n+l-i)/21/T( (n+hn+l-i) /2 ) 


(14) 
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Korin  [17]  computed  percentage  points  by  using  Box's  asymptotic  ex- 
pression up  to  terms  of  order  n Davis  and  Field  [8]  computed  the 

percentage  points  by  using  the  Cornish-Fisher  type  inversion  of  Box's 
asymptotic  expression.  In  this  paper,  we  approximated  the  distribution 
of  V^b,  with  a suitable  Pearson  type  distribution  where  b is  a suitably 
chosen  integer.  The  type  of  distribution  that  is  fitted  is  the  Pearson's 
Type  I distribution.  In  Table  VI,  we  compared  some  of  the  percentage 
points  obtained  by  us  with  the  corresponding  values  obtained  by  Davis 

and  Field  [8]  and  Korin  [17],  where  n.  = n . The  percentage  points 

1 o 

c4  obtained  by  us  using  the  Pearson  type  approximation  are  given  under 
the  column  L-OK  whereas  the  corresponding  values  obtained  by  Korin,  and 
Davis-Field  are  given  under  the  columns  Korin  and  D-F,  respectively. 

Table  VI  indicates  that  the  Pearson  type  approximation  is  satisfactory 
for  practical  purposes.  In  this  table,  a is  defined  by  the  equation 

P [ i c4IH4  ] =(!-“>  <15) 

When  n is  large,  Bishop  [18]  investigated  the  accuracy  of  approximating 
o 

2 

the  distribution  of  the  -th  power  of  the  likelihood  ratio  criterion  with 

N 

the  Beta  distribution  for  a few  cases,  where  n * qn  . 

o 

Using  the  Pearson  type  approximation,  we  computed  the  values  of  c4 
when  nQ  = (p+1) (1) 20(5) 30,  p * 2(1)6,  q = 2(1)10  and  a = 0.01,  0.05, 

0.10.  These  percentage  points  are  given  in  Table  VII. 
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There  is  a typographical  error  in  Korin’s  table.  The  correct  value  seems  to  be  39.65 


TABLE  VII 


PERCENTAGE  POINTS  OF  THE  DISTRIBUTION 
OF  -2  log  V4 


a = o.io  p = 2 


n 

o 

1 

to 

q*3 

q«4 

q=5 

q=6 

>Q 

II 

"■j 

q«8 

q=9 

3 

9.75 

15.80 

21.29 

26.52 

3l.o0 

3b. 58 

41.49 

«*6.32 

4 

0.56 

14.07 

19.08 

23.87 

28.52 

33.17 

37.55 

41.97 

5 

7.97 

13.22 

17.99 

22.55 

26.  J7 

31.31 

35.50 

39.80 

6 

7.63 

12.  70 

17.33 

21.75 

26.05 

30.25 

3*  • 39 

38.49 

7 

7.39 

12.36 

16.  09 

21.22 

25. <*2 

29.54 

33.60 

37.61 

0 

7.23 

12*11 

16.57 

20.84 

24.98 

29.04 

33.  03 

36.90 

9 

7.10 

11*93 

16.34 

20.55 

24.55 

28.  b6 

32.61 

36.51 

10 

7.01 

11.70 

16.16 

20.33 

2** . J9 

28.36 

32.28 

36.14 

li 

6.93 

11.67 

16.01 

20.15 

24,18 

28.12 

32.01 

35.05 

12 

6.87 

11.50 

15.00 

20.01 

24.01 

27.93 

31.79 

35.61 

13 

6.02 

11.50 

15.79 

19.89 

23.07 

27.77 

31.62 

35.41 

1 4 

6.77 

11.43 

15.70 

19.78 

23.75 

27.63 

31.46 

35.24 

15 

6*74 

11.37 

15.63 

19.69 

23.65 

27.52 

31.33 

35.  09 

16 

6.71 

11.  33 

15.57 

19.62 

23.56 

27.42 

31.21 

34.97 

17 

6.68 

11*28 

15.51 

19.55 

23. 4 0 

27.32 

31.11 

34.06 

10 

6.65 

11.25 

15*46 

19.49 

23..1 

27.25 

31.03 

34. 7o 

19 

6.63 

11.21 

15*42 

19.44 

2?.35 

27.18 

30.95 

34.60 

20 

6.61 

11.10 

15.30 

19.40 

23.30 

27.12 

30.88 

34.60 

25 

6.53 

11*07 

15.24 

19.22 

23.09 

26.  06 

30*61 

34.30 

30 

6.40 

11*00 

15.14 

19.10 

22.96 

26.73 

30.45 

34. 12 

a = 0. 

10  p 

= 3 

n 

o 

qss2 

q-3 

q*4 

q-s 

q=6 

q«7 

q-8 

q=9 

4 

10.07 

30.76 

41.70 

52.41 

62*90 

73.03 

03*14 

93.  lb 

5 

16.19 

26.87 

36.79 

46.36 

55.73 

64.95 

74.05 

83.  09 

6 

14.03 

24.07 

34.22 

43.24 

52.06 

60.76 

b9  » 35 

^7.86 

7 

14.01 

23.66 

32*64 

41.33 

45.91 

58.10 

66  • *4 

/ 4 ,64 

0 

13.46 

22.03 

31.56 

40.03 

40.29 

5b.  44 

64 » *0 

72.45 

9 

13*06 

22.24 

30.51 

39.09 

47.19 

55.17 

63  • 0o 

70.00 

10 

12.77 

21.  79 

30.23 

38.38 

4b.  3b 

54.21 

61.90 

69.68 

11 

12.53 

21.44 

29.77 

37.82 

45.70 

53.  4b 

61.13 

66.74 

12 

12.35 

21.16 

29.40 

37.37 

' 5.17 

52.06 

b0  • 4j 

67.98 

13 

12.19 

20.  93 

29.10 

37.00 

44 . 7h 

52.  3b 

59.0* 

67.35 

14 

12.07 

20.73 

28.95 

36.69 

44.3b 

51.94 

59.42 

66.84 

15 

11.96 

20.57 

28.63 

36.43 

44.36 

51.59 

59.03 

66.39 

lb 

11.07 

20.43 

28.45 

36.21 

43.90 

51.20 

50. o7 

6b.  0< 

17 

11*70 

20.  30 

20.29 

36.01 

43,57 

51.02 

50.30 

65  .b7 

10 

11.71 

20.20 

20.15 

35.83 

43.  37 

50.79 

50.11 

65.30 

19 

11*65 

20.10 

20.02 

35.68 

43.19 

50.58 

57.08 

65.12 

20 

11*60 

20.02 

27.91 

35.55 

*>3.02 

50. 39 

57  *60 

64.90 

25 

11*39 

19.70 

27.50 

35.05 

42.43 

49.71 

56.91 

64.04 

30 

11.26 

19.50 

27.23 

34.72 

42.05 

49.27 

5b. 42 

63.49 

q=10 


51.12 
4b  .36 

43.97 
42.53 
41.57 
4 0.60 

40.37 

39.97 
39.65 

39.38 
39.16 

38.98 
38.82 
30.00 
30.56 
30.45 
30*36 
36.20 
37.95 
37.74 


q=10 


103*12 

92.05 

86.30 

02.76 
80.  37 
78.64 

77.31 
7b. 2& 
75.45 

74.77 
7-.  19 
73.70 
73.29 
7 2.92 
72.59 

72.32 
72. G7 
71.12 

70.32 


22 


TABLE  VII  (Continued) 


a = 0.10  p = 4 


n 

o 

q=2 

q*3 

q=4 

q=5 

q=6 

q=7 

q=6 

q=9 

q=10 

5 

30.70 

50.45 

68.96 

86.91 

104.53 

121.93 

139.15 

156.25 

173.24 

o 

26.16 

43.85 

60.44 

76.56 

92.38 

107.98 

123.4o 

136.81 

154. 09 

7 

23.79 

40.  34 

55.91 

71.02 

85.37 

100.51 

115.03 

129.45 

143.75 

3 

22.33 

38.15 

53.  06 

67.55 

81.76 

95.  81 

109.73 

123.52 

137.25 

9 

21*33 

36.65 

51.  11 

65,15 

78.  J4 

92.55 

106. 0o 

119.46 

132.76 

10 

20.60 

35.56 

49,67 

63.39 

76.37 

90,17 

103.37 

116.45 

129.48 

11 

20.06 

34.72 

46*58 

62.04 

75.28 

88.36 

101. 30 

114.15 

126.93 

12 

19.62 

34.06 

47.71 

60.96 

74,02 

86*90 

99.66 

112.34 

124.95 

13 

19.27 

33.53 

47.00 

60.12 

73.00 

85.73 

96.36 

110. 66 

123.31 

14 

18.93 

33k  08 

46.42 

39.40 

72,17 

84.  7b 

97.26 

109.65 

121.98 

15 

18. 7h 

32.71 

45.94 

56.80 

71.45 

83.95 

96.32 

108.61 

120.85 

16 

18.53 

32.40 

45.52 

58.30 

70.34 

83.25 

95.55 

107*73 

119.87 

17 

16.36 

32.12 

45,16 

57,65 

70.33 

82.65 

94.86 

106.98 

119.03 

18 

18.21 

31.89 

44,85 

57.47 

65.97 

82.13 

94.26 

106.32 

116.30 

19 

18.07 

31.68 

*4.58 

57.13 

69.47 

81.66 

93.75 

105.73 

117.67 

20 

17,96 

31.49 

44.33 

56.83 

69.12 

81.  25 

93.28 

105.22 

117.10 

25 

17.51 

30.81 

43.  *.3 

55.72 

67.30 

79.75 

91.57 

103.33 

115.00 

30 

17.24 

30.37 

42.86 

55.02 

66.37 

78.79 

90.49 

102. 11 

113,67 

a 

= 0.10  p 

= 5 

n 

O 

q«2 

q=3 

q=4 

q-5 

q=6 

q=7 

q-8 

q=9 

q=10 

6 

45.25 

74.90 

102.86 

130.09 

156,89 

183.39 

209.65 

235.79 

261.72 

7 

38.56 

65.12 

90.21 

114.66 

138. 7 0 

l<62, 5 0 

186.15 

209.58 

232*90 

8 

34,95 

59.77 

83.27 

106.16 

128.  70 

150. 99 

173.1? 

195. 06 

216.94 

9 

32.69 

56.37 

78.81 

100.70 

122.23 

143.54 

164.69 

185.70 

206.61 

10 

31.12 

54.00 

75.70 

96.88 

117. 72 

138. 36 

! 56  * 82 

179.16 

199.37 

11 

29.97 

52.24 

73.40 

94.03 

114,37 

134. 47 

154.42 

174.28 

194.00 

12 

29.08 

50.38 

71.62 

91.85 

111.78 

131.  47 

151.03 

170.52 

189.84 

13 

28.39 

49.31 

70.20 

90.09 

109.70 

129.10 

148.33 

167.49 

186.93 

14 

27.82 

48.93 

69.04 

88.66 

108.01 

127. 14 

146.14 

165. 04 

183.78 

15 

27,35 

48.21 

68.07 

87,48 

106,60 

125.51 

1h4 ■ 31 

162.94 

181.54 

16 

26.95 

47.59 

67,26 

86.47 

105.42 

124.14 

142.75 

161.24 

179.62 

17 

26.62 

47.07 

66.57 

85.62 

104.40 

122.98 

141.42 

159.75 

177.97 

18 

26.32 

46.61 

65.  96 

84.87 

103.30 

121.95 

140.27 

158.45 

17b.  56 

19 

26.07 

46.21 

65.44 

84.22 

102.74 

121.06 

139,23 

157.32 

175. 30 

20 

25.85 

45.  35 

64.97 

83,65 

102.05 

120.28 

138.36 

156.32 

174. 20 

25 

25.03 

44,59 

63*28 

81,56 

95,57 

117.  41 

135.10 

152.67 

170, 19 

30 

24.  h> 

43.78 

62.22 

80.23 

98.01 

115. 60 

133.04 

150.39 

167,69 

23 


TABLE  VII  (Continued) 


a * O.IO  p = 6 


n 

0 

q=2 

q=*3 

q=4 

q=5 

q=*6 

q-7 

q=6 

q=9 

q=10 

7 

62  + 54 

104.12 

143.53 

181.97 

215.97 

257.40 

294.64 

331.76 

368.56 

8 

53. *2 

90.76 

126.20 

160.82 

194.34 

228.72 

262.28 

295.59 

328.81 

9 

4tt.n0 

63.29 

110.44 

146.84 

180.79 

212.41 

243.33 

275.04 

306*13 

10 

45.18 

76.42 

110,07 

141.01 

171.49 

201.73 

231.72 

261.56 

291.,  26 

11 

*2.92 

74.99 

105.56 

135.44 

164.31 

194. 13 

223.14 

251.93 

260.69 

u 

41.23 

72.42 

102.16 

131.26 

159.38 

188.42 

216.65 

244.  70 

272,69 

13 

39.94 

70.43 

99.5,4 

128.00 

156.11 

193.93 

211.57 

239.10 

266.46 

14 

38.91 

68.83 

97.42 

125.39 

153.01 

180,35 

207.54 

234.56 

261.4b 

15 

3b. 06 

67.51 

95.67 

123.24 

150. 46 

177. 39 

204.16 

230.76 

257.33 

16 

37.  $6 

66.  -.2 

94.23 

121.46 

148. 31 

174.95 

201.33 

227.67 

253.85 

17 

36.76 

65.49 

93.00 

119.93 

146.51 

172.84 

199.02 

225.05 

250.97 

Id 

36.26 

64.69 

91.93 

116.60 

144.95 

171.03 

196,97 

222.74 

246. *4 

19 

35.82 

64.00 

91.01 

117.47 

143.59 

169.49 

195.17 

220.73 

246.20 

20 

35.43 

63.39 

90.21 

116.47 

142.40 

166.11 

193.61 

219.00 

244.29 

25 

34.06 

61.22 

87.30 

112.66 

138.11 

163.13 

187.96 

2 1C • 70 

237*33 

30 

33.21 

59.97 

85.49 

110.63 

135.45 

160.06 

184.51 

208.78 

233.04 

a « 0.05  p « 2 


n 

0 

q-2 

q=3 

q=4 

q«5 

q«6 

q-7 

q=8 

q*9 

q=10 

3 

12.18 

18.70 

24.55 

30.09 

J5..5 

40.66 

45.81 

50.87 

55.86 

4 

10.70 

16.65 

22.00 

27.07 

31.  37 

36.75 

41  * 45 

4b. 07 

50.  b4 

5 

9.97 

15.63 

20.7  3 

25.57 

30.23 

34.79 

39.26 

43.67 

48.02 

6 

9.53 

15.02 

19.97 

24.66 

29.19 

33.61 

37.95 

42.22 

46.45 

7 

9.24 

14.62 

19.46 

24.05 

28.49 

32.83 

37. 08 

41.2b 

45.40 

3 

9.04 

14.33 

19.10 

23.62 

22.99 

32.26 

36.44 

40.57 

44.64 

9 

8.88 

14.11 

16.83 

23.30 

27.62 

31.84 

35.98 

40.05 

44.08 

10 

8.76 

13.94 

16.61 

23.05 

27.  J3 

31.51 

35.61 

39.65 

43.64 

11 

8.67 

13.  31 

18.44 

22.65 

27.10 

31.25 

35.32 

39.33 

43.29 

12 

8.59 

13.70 

16.30 

22.66 

26.90 

31.03 

35.  Oi 

39.07 

43.  OC 

13 

8.52 

13.60 

18.  19 

22.54 

26.75 

30.  65 

34.67 

36.64 

42.76 

14 

8,47 

13.53 

18.10 

22.42 

26, g 1 

30.70 

34.71 

36.  cb 

42,56 

15 

8.42 

13.46 

18.01 

22.33 

26.50 

30.57 

3n  ,57 

36.50 

42.36 

16 

6,36 

13.40 

17,94 

22.24 

26. n 0 

30.45 

34.43 

38.36 

42.23 

17 

8.35 

13.35 

17.87 

22.17 

26.31 

30.  35 

34.32 

36.24 

42.  10 

16 

6.32 

13.30 

17.82 

22.10 

26.23 

30.  27 

34.2  3 

36.13 

41.99 

19 

6.26 

13.26 

17.77 

22.04 

26.16 

30.19 

34.14 

36.04 

41.86 

20 

6.26 

13.23 

17.72 

21.96 

26.10 

30.12 

74,07 

37.95 

41.79 

25 

6.17 

13.10 

17.55 

21.79 

25.87 

29.66 

33.76 

37.63 

41.44 

30 

8.11 

13.01 

17.44 

21.65 

25.72 

29.69 

33.59 

37.42 

41.21 

24 


TABLE  VII  (Continued) 


a = 0.05  p - 3i 


n 

o 

q«2 

q=3 

q=4 

q*5 

q=6 

q=7 

q=8 

q=9 

q=10 

4 

22.41 

35,00 

46,58 

57.68 

68.50 

79.11 

89.60 

99.94 

110.21 

5 

19.19 

30.52 

40.95 

50.95 

60.59 

70.26 

79.69 

89.03 

98.27 

6 

17,57 

28.24 

38.06 

47.49 

56.o7 

65.69 

74.58 

83.39 

92.09 

7 

16.59 

26.34 

36.29 

45.37 

54.20 

62.89 

71.44 

79.90 

88.30 

8 

15.93 

25.90 

35.10 

43.93 

52,5*. 

60.99 

69.32 

77.57 

85.73 

9 

15.46 

25.22 

34.24 

42.90 

51.33 

59.62 

67.78 

75.86 

83.87 

19 

15.11 

24.71 

33.59 

42.11 

50,42 

58.57 

66.62 

74.58 

82.46 

11 

14.83 

24.31 

33.08 

41.50 

49.71 

57.76 

65.71 

73.57 

81,36 

12 

14.61 

23.99 

32.67 

41.00 

49.13 

57.11 

64.97 

72.75 

80.45 

13 

14.43 

23.73 

32.33 

40.60 

48.65 

56.56 

64.36 

72.09 

79.72 

14 

14.28 

23.50 

32.05 

40.26 

48.26 

56.11 

63.86 

71.53 

79.11 

15 

14.15 

23.32 

31.81 

39.97 

47.92 

55.73 

63.43 

71.05 

78.60 

16 

14.04 

23.16 

31.60 

39.72 

47.63 

55.40 

63.06 

70.64 

78.14 

17 

13.94 

23.02 

31.  h3 

39.50 

47.38 

55,11 

62.73 

70.27 

77.76 

18 

13,86 

22.39 

31.26 

39.31 

47.16 

54,86 

62.49 

69.97 

77.41 

19 

13.79 

22.78 

31.13 

39. *5 

46.96 

54.64 

62.21 

69.69 

77.11 

23 

13.72 

22.69 

31.01 

39.00 

46.79 

54.44 

61.98 

69. h5 

7o.  64 

25 

13.48 

22.33 

30.55 

38.44 

46*15 

53,70 

61.16 

68.54 

75.84 

30 

13.32 

22.10 

30.25 

38,09 

45.73 

53.22 

60.62 

67.94 

75.18 

a = 

0.05  p 

= 4 

n 

O 

q=2 

q=3 

q-4 

q-5 

q=6 

q-7 

q=8 

q=9 

q=10 

5 

35.39 

56.10 

75.36 

93.97 

112.17 

130.11 

147.81 

165.39 

192.60 

6 

30.06 

48.62 

65.90 

82.60 

96.93 

115.03 

130.94 

146.69 

162.34 

7 

27.31 

44.69 

60.89 

76.56 

91.86 

106.96 

121.90 

136.71 

151.39 

8 

25.61 

42.24 

57.77 

72.77 

87 « >.6 

101.94 

116.23 

130.43 

144.50 

9 

24.45 

40.57 

55.62 

70.17 

64.42 

96.46 

112.32 

126.08 

139.74 

10 

23.62 

39.34 

54.04 

68.26 

82.19 

95.90 

109.46 

122.91 

136.24 

11 

22.98 

38*41 

52.84 

66.81 

80.48 

93.95 

107.27 

120.46 

133.57 

12 

22.48 

37.67 

51.90 

65,66 

79.14 

92.41 

105.54 

116.55 

131. *5 

13 

22.08 

37.08 

51.13 

64.73 

78.04 

91.15 

104.12 

116.98 

129.74 

14 

21.75 

36.59 

50.50 

63.95 

77.13 

90. 12 

102.97 

115.69 

128.32 

15 

21.47 

36.17 

49.97 

63.30 

76.37 

89.26 

101. 99 

114.69 

127.14 

16 

21.24 

35.32 

49.51 

62.76 

75.73 

88.51 

101.14 

113.67 

12o*10 

17 

21.03 

35.52 

49.12 

62.28 

75.16 

87.87 

100.42 

112.67 

126.22 

18 

20.86 

35.26 

48.78 

61*66 

74.68 

87.  31 

99,80 

112.17 

124.46 

19 

20.70 

35.02 

48.47 

61.50 

74.25 

86.82 

99.25 

111.56 

123.79 

20 

20.56 

34.82 

48.21 

61.17 

73.87 

86.38 

98.75 

111.02 

123.16 

25 

20.06 

34.06 

47.23 

59.98 

72.47 

84c  78 

9b.  95 

109.01 

120,99 

30 

19.74 

33.59 

46.61 

59.21 

71.56 

83.74 

95.79 

107.71 

119.57 

25 


TABLE  VII  (Continued) 


a = 0.05  p * 5 


n 

o 

q«2 

q-3 

q*4 

q-5 

q=6 

q®? 

q=8 

q=9 

q=10 

6 

51.11 

01.99 

110.92 

130.90 

166.54 

193.71 

220.66 

247.37 

273.  00 

7 

43.40 

71.06 

97*03 

122.22 

146.95 

171. 34 

195.43 

219.47 

243.30 

6 

39.29 

65.15 

09.45 

113.03 

136.16 

159.04 

101.65 

204.14 

226.48 

9 

36.71 

61.39 

04*62 

107.17 

129.30 

151.17 

172.80 

194.27 

215.64 

10 

34.93 

50.70 

01.25 

103.06 

124.40 

145. 64 

166.5b 

107.37 

208.02 

11 

33.62 

56.05 

70.75 

100.02 

120.92 

141. 54 

161*93 

102.24 

202.37 

12 

32.62 

55.37 

7b. 03 

97.60 

110.15 

138.30 

158.38 

178.23 

198.03 

13 

31.03 

54.19 

75.30 

95.02 

115.36 

135.86 

15s .54 

175.10 

194.51 

m 

31.19 

53.23 

74.05 

94.29 

114.16 

133.80 

153.21 

172.49 

1 91  * t>8 

i5 

30.66 

52.44 

73.01 

93.02 

112.56 

132. 07 

151.29 

170.36 

189.30 

16 

30.22 

51.76 

72.14 

91.94 

111.41 

130. 61 

149. 6o 

1 op . 53 

167.32 

17 

29.03 

51.19 

71.39 

91.03 

110.34 

129. 38 

148.25 

166.99 

185.61 

10 

29.51 

50.69 

70.74 

90.23 

109.39 

128.29 

147.03 

165.65 

184.10 

19 

29.22 

50.26 

70.17 

09.54 

100.57 

127. 36 

145.97 

164.45 

182.81 

20 

20.97 

49.00 

69.67 

00.93 

107.35 

126.52 

145.02 

163.38 

181.65 

25 

20.05 

40.40 

67.06 

06.70 

105.21 

123.51 

141.62 

159.60 

177.49 

30 

27.40 

47.61 

66.71 

05.29 

103.56 

121.60 

139,47 

157.22 

174. b7 

a 

* 0.05  p 

= 6 

n 

o 

q=2 

q=3 

q=4 

q=5 

q=6 

q=7 

q=8 

q=9 

q=10 

7 

69.62 

112.69 

153.28 

192.75 

23l.o0 

269.98 

308.03 

345. 81 

333.45 

6 

59.24 

97.92 

134.42 

169.94 

204. 36 

239.  37 

273.52 

307.01 

341.44 

9 

53.5b 

89.72 

123.87 

157.11 

189.33 

222.14 

254.18 

286.00 

317.69 

10 

49.95 

84.42 

117.03 

i4S.7r 

179. 37 

210.85 

241.48 

271.84 

30?.  n 

11 

47.43 

80.69 

112. 18 

142.85 

173.  10 

202. 35 

232.46 

261.80 

291.03 

12 

45.56 

77.90 

108.56 

138.40 

167.79 

196.84 

225.65 

254.25 

282.76 

13 

44.11 

75.74 

105.72 

134.94 

163.72 

192.15 

220.37 

248.39 

276.24 

14 

42.96 

74.00 

103.46 

132,19 

160.44 

188. 39 

216.09 

243. b7 

271.04 

15 

42.03 

72.59 

101.61 

129.90 

157.76 

185.31 

212.63 

239.72 

266.75 

16 

41.24 

71.41 

100.07 

128,01 

155.52 

182.73 

209.70 

236.48 

263.15 

17 

40.50 

70.40 

98.76 

126.38 

153,60 

180.51 

207.23 

233.73 

260. p7 

18 

40.02 

69.54 

97.63 

125.01 

151. 36 

178.63 

205.09 

231.36 

257.48 

19 

39.53 

60.30 

96.64 

123.80 

150.55 

176.99 

203.21 

229.26 

.55.19 

20 

39.11 

68.13 

95.77 

122.73 

149.29 

175.53 

201. oO 

227.45 

253.19 

25 

37.58 

65.79 

92.69 

118.93 

144.78 

170.35 

195.69 

220.92 

246.00 

30 

36.63 

64.33 

90.77 

116;57 

141. 37 

167.13 

192.04 

?lo.82 

241.49 

26 


TABLE  VII  (Continued) 


i 

i 

! 


a = o.oi 


V 

f 


n q-2 

o 

3 17.67 

4 16.53 

5 14.40 

6 13.84 

7 13.42 

8 13.12 

9 12.90 

10  12.72 

11  12.58 

12  12.47 

13  12.37 

14  12.30 

15  12.23 

16  12.17 

17  12.12 

18  12< 07 

19  12.03 

20  11.99 

25  11.86 

30  11.77 


q=*3  q-4 


24.99  31.49 
22.25  28.20 

20.38  26.56 
20.07  25.58 
19.52  24.93 
19.13  24.46 
13*.  34  24.11 

13.ol  23.84 

18.43  23.62 
18.29  23.44 

18.16  23.30 
18.06  23.17 
17.97  23.07 
1' • 89  22.97 
17.83  22.89 
17.77  22.82 
17.71  22.75 
17.67  22.69 
17.49  22.48 
17.  37  22.  34 


q=5 

q=6 

37.60 

43.45 

33.80 

39.15 

31.90 

37.01 

30.76 

35.72 

30.00 

34.16 

29.46 

34.25 

29.05 

33.80 

28.74 

33.44 

26.49 

33.15 

26.28 

32.92 

26.11 

32.72 

27.97 

32.56 

27.64 

32.42 

27.73 

32.29 

27.64 

32.18 

27,56 

32.39 

27.46 

32.01 

27.42 

31.93 

27.16 

31. i? 

27.00 

31. *7 

q=7  q=8 


49.14  54.70 
44.36  49.44 
41.97  4b • 82 
40.54  45.24 

39.59  44.20 
38.90  43.4-* 

38.39  42.83 
37.99  42 . 4-. 
37.68  42.09 

37.41  41.80 
37.20  4i.57 

37.01  41.36 
36.85  41.19 
36.72  41.04 

36.60  40.91 
36.  49  40.79 

36.40  40.69 
36.  32  40.60 
36.00  40.25 
35.80  40.02 


q=9  q=10 


60.16  65.52 

54.43  59.35 
51.57  56.26 
49.86  54.40 
48.71  53.16 

47.89  52.28 
47.28  51. bl 
46.60  51.10 

46.42  50.69 
46.11  50.35 
45.85  50.07 
45.63  49.83 

45.44  49.62 
45.26  49.45 
45.13  49.29 

45.01  49.16 

44.89  49.03 
44.79  48.93 

44.42  46.52 

44.17  46.25 


a * 0.01  p = 3 


n 

O 

q-2 

q=3 

q-4 

q=»5 

q=6 

q-7 

q=8 

q*9 

q=10 

4 

30.17 

44.0  0 

56.58 

66.55 

80.14 

91.49 

102.60 

113.57 

124.41 

5 

25,75 

36.23 

49.56 

60.36 

70.  *0 

81.01 

91.03 

100.91 

110.66 

6 

23.53 

35.30 

46.00 

56.17 

66.02 

75.65 

85.10 

94.42 

103.61 

7 

22.20 

33.53 

43.83 

53.62 

63.11 

72.36 

61.47 

90.45 

99.30 

6 

21.31 

32.34 

42.37 

51.91 

61.14 

70.17 

79.03 

87.76 

96.36 

9 

20.67 

31.46 

41.32 

50.67 

59.73 

69.58 

77.25 

65.62 

94.27 

10 

20.19 

30.84 

40.52 

49.73 

56.65 

67.37 

75.91 

64.34 

92.69 

11 

19.62 

30.  33 

39.90 

49.01 

57.81 

66.42 

74.66 

63.20 

91.43 

12 

19.52 

29.93 

39.40 

46.42 

57.14 

6^.66 

74.03 

62.26 

90.  **2 

13 

19.26 

29.50 

39.00 

47.93 

5t.?6 

65.  03 

73.33 

81.61 

6 1.59 

1 4 

19. 0’ 

29.32 

36.65 

47.63 

56.13 

64.51 

72.75 

60.67 

98,90 

15 

16.90 

29.0a 

36.36 

47.16 

55. ’3 

64.07 

72.27 

60.34 

68.31 

16 

16.75 

26.86 

39.11 

46.69 

55,39 

63.69 

71.83 

79.6c 

87.75 

17 

16.62 

26.71 

37.  89 

46.63 

55.10 

63.36 

71.47 

79.47 

67.36 

16 

18.51 

28,55 

37.70 

46.41 

54.34 

63.07 

71.14 

79.12 

66.97 

19 

16. Hi 

28.42 

37.54 

46.21 

94.61 

62.60 

70,87 

78.60 

66.64 

20 

16.32 

26.30 

37,39 

46.04 

54. *1 

62.56 

70.60 

76.62 

6c.  33 

25 

17.99 

27.85 

36.63 

45.36 

53.^5 

61.74 

69.66 

77.48 

65.20 

30 

17.79 

27.56 

36.47 

44.96 

53.16 

61.18 

69.05 

76.61 

64.46 

27 


TABLE  VII  (Continued) 


a = O.Ol  p = 4 


no 

q=2 

q-3 

q=4 

q=s 

q*6 

q=7 

q=8 

q=9 

q=10 

5 

45.49 

67.39 

99.52 

109.33 

127.60 

146.51 

165.13 

183.54 

201.62 

6 

39.41 

59.52 

77.05 

94.91 

112.11 

129. 07 

145.79 

162.32 

176.70 

7 

34.90 

53.66 

71.05 

87.72 

103.95 

119.66 

135.53 

151.06 

166.42 

8 

32.59 

50. 06 

67.33 
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95.36 
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153.43 
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142.39 
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27.60 
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59.74 

73.09 

87.05 
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114.26 

127.60 

140.81 
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59.12 

72.34 

86.19 

99.76 

113.16 
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139.49 
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71.71 
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125.35 

138.36 

17 
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42.50 

57.13 

71.17 

84.62 

96.23 
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137.41 
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56.72 

70.68 

64.26 

97,59 
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136.56 

19 
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56.39 

70.26 

83.79 

97.04 
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135.61 

20 
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41.65 

56.06 

69.89 

63.36 

96.55 

109.56 

122.42 

135.13 

25 
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40.74 

54.92 

68.52 

61,77 

94.75 

107.55 

120.19 

132.73 

30 
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40.17 

54.19 

67.64 
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93.60 

106.26 

119.78 
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a * 0.01  p * 5 
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10 
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126.32 
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169.66 
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14 
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61.95 

64.10 
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146.66 
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167.15 
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15 
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162.66 
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17 
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161.67 
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100.66 
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198.86 
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100.10 
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139.76 

159.14 
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20 
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79.09 

99.41 

119.29 

136.62 

156.11 

177.22 

196.16 

25 

34,33 

56.36 

77.00 

96.91 

116.36 

135.49 

154.36 

173.09 

191.64 

30 

33.62 

55.37 

75.70 

95.34 

114.52 

133.40 

152.02 

170.46 

166.62 
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a = o.oi  p = 6 
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307.34 

340.08 

10 

59.81 
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196.59 

228.  79 

260.53 
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11 
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92.15 

125.34 

157.47 

188.32 

219.93 

250.64 

280.98 

311.19 

12 

54.43 

fed. 92 

121.24 

152.51 

193.14 

213.36 
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272.79 

302.24 

13 

52.68 

8o . 4 1 

118.04 

148.67 

178.64 

208.  24 

237*44 

266.46 

295.22 

14 

51.29 

84.41 

115.48 

145.56 

175.06 

204. 06 

232.85 

261.33 

289.58 

15 
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82.77 

113.38 

143*04 
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284.98 

16 

49.22 

81.  -*1 

111.65 

140.92 

169.61 
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253.59 
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17 
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19 

*♦7,15 

78.40 
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272.61 

20 
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190. 07 
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25 
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103. 36 
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157.37 

184.  44 

21C.76 

236.82 

262.71 

30 

43.67 

73.30 

101.21 

128.26 

154.30 

180.93 

206.79 

2 32.45 

257.88 
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SECTION  VII 


LIKELIHOOD  RATIO  TEST  STATISTIC  FOR  MULTIPLE 
HOMOGENEITY  OF  THE  COVARIANCE  MATRICES 

In  this  section,  we  discuss  Pearson  type  approximations  to  the  dis- 
tribution of  the  likelihood  ratio  statistic  for  testing  the  hypothesis 
when  p^  = p(i  = l,...,q)  and  where  H5  was  defined  in  Section  II.  This 
hypothesis  is  of  interest  in  studying  certain  linear  structures  on  the 
covariance  matrices  (see  Krishnaiah  and  Lee  [19]).  The  likelihood  ratio 
statistic  for  testing  H^_  is 
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where  ni  and  were  defined  in  the  preceding  section,  and 
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The  h-th  moment  of  is  given  by 
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As  in  the  preceding  section,  the  distribution  of  is  approximated 

with  the  Pearson's  Type  I distribution  where  b is  a properly  chosen 

integer.  In  Tables  VIII  and  IX,  we  have  n^  = n^  and  q - kd. 

In  Table  VIII,  a^  is  the  value  of  a obtained  if  we  use  Pearson 

type  approximation  whereas  cl^  is  the  value  of  a obtained  by  using  Box's 

-13 

asymptotic  series  up  to  order  n , where 


(l-a) 


(18) 


and 


V*  = -2  log  V5. 


Table  VIII  indicates  that  the  accuracy  of  the  Pearson  type  approximation 
is  sufficient  for  practical  purposes.  Using  the  above  approximation,  ..e 
computed  tables  for  the  values  of  where  a = 0.05,  0.01,  k ■ 2,  3, 

M = 1(1)20(5)30,  M = nQ  - p,  p * 1,2, 3, 4.  These  values  are  given  in 
Table  IX. 
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APPPTOX 


Let  W (0  < W < 1)  be  a random  variable  whose  moments  are  given  by 


e(w‘j  = c[-fr  y//txf‘)4 

C H Xl 


if  r t 


tt  r f;J 


0 a u 

where  C is  chosen  such  that  E (Wj  = 1 and  J X.  = £ Y . . Box  (1949)  gave  an 

i=l  1 j=l  J 

asymptotic  expression  for  the  distribution  of  -2p  log  W,  where  p (0  ^ p ^ 1)  is 
an  arbitrary  constant. 

In  this  appendix  we  give  asymptotic  expression  for  the  distribution  of 
-2P  log  W with  remaining  term  of  order  where  p is  chosen  such  that  = 0, 
wr  is  given  by 


•'  Ik  (FxTT~  % “( ftrl 

= Cl  “ p)^,  z.  - (1  - p)Y^  and  Br(h)  is  the  Bernoulli  polynomial  of  degree 
r and  order  unity. 

Let  gv(z)  be  the  probability  density  function  (p.d.f.)  of  the  x distribution 
with  v degrees  of  freedom  and  N be  the  sample  size.  'Hie  p.d.*\  of  -2p  log  W is 
given  by 


im  = £.  SrC})  4 OL*"r> 
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where 
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